function [m, S, deg2rad, R0, g0, Vc, rho0, hs, sigmamin, sigmamax, dQmax, CQ, qmax, nmax, kQ] = getConstant()
    m = 104305;
    S = 391.22;
    deg2rad = pi/180;
    R0 = 6378137;
    g0 = 9.81;
    Vc = sqrt(R0*g0);
    rho0 = 1.2256;
    hs = 7254.24;
    sigmamin = 0;
    sigmamax = 100*deg2rad;
    dQmax = 15e5;
    CQ = 1e-8;
    qmax = 18e3;
    nmax = 2.5;
    kQ = 9.4369e-5;
end

